On stability of odd-frequency superconducting state 
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Odd-frequency pairing mechanism of superconductivity has been investigated for several decades. 
Nevertheless, its properties, including the thermodynamic stability, have remained unclear. In 
particular, it has been argued that the odd-frequency state is thermodynamically unstable, has an 
unphysical (anti-) Meissner effect, and thus can not exist as a homogeneous equilibrium phase. We 
argue that this conclusion is incorrect because it implicitly relies on the inappropriate assumption 
that the odd-frequency superconductor can be described by an effective Hamiltonian that breaks 
the particle conservation symmetry. We demonstrate that the odd-frequency state can be properly 
described within the functional integral approach using non-local-in-time effective action. Within the 
saddle point approximation, we find that this phase is thermodynamically stable, exhibits ordinary 
Meissner effect, and therefore can be realized as an equilibrium homogenous state of matter. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.20.Rp 
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Odd-frequency pairing mechanism was first introduced 
by Berezinski:^ in an attempt to describe the A-phase of 
superfluid "^He. Although, it was latter found that in 
the case of superfluid '^He the pairing is odd is space 
(p-wave) rather than in time, the question weather the 
same-spin fermions can develop anomalous correlations 
that are odd in frequency ( "p-wave" along the time direc- 
tion) remained open. The interest to this pairing mech- 
anism was revived almost two decades later in the con- 
text of high-temperature superconductors by Balatsky 
and Abrahams,^ who generalized it to the case of oppo- 
site spin pairing. At the same time it has been argued 
that the odd-frequency superconducting order parame- 
ter does not correspond to a thermodynamically stable 
phase, and thus may not be realized as a homogenous 
equilibrium state of matteri^i^ On the other hand, the 
authors of Ref. have shown that path integral formula- 
tion of the problem can lead to the opposite conclusion. 
No resolution of this contradiction have been proposed 
so far. In this report we demonstrate that odd-frequency 
superconductor is thermodynamically stable and exhibits 
ordinary Meissner effect. 

The argument on instability of the odd-frequency su- 
perconductor can be summarized as follows'^''*. Thermo- 
dynamic stability is determined by the change of free en- 
ergy in the ordered phase, relative to the disordered one. 
Near the second-order transition point it can be written 



AO 1 A(a.,q)A+(^,q) 



(1) 



with gap-functions A(w, q) and A+(ti),q) related 
to the anomalous (Matsubara) Green's functions 
F(w,q) = /dre^'^^(T^Cq(T)c_q(0)) and i^+(w,q) = 

/ (iTe*'^'^(TT-cLq(T)c^(0)) through a self-consistency re- 
lation 



D{u;-u;',q-q')F{u;',q'), (2) 



and an identical relation for and A+ . Here D{uj,q) 
is irreducible interaction between quasiparticles with en- 
ergies (we assume that D is real and even in both 
uj and q), and /3 is the inverse temperature. Here and 
in the following we consider spinless fermions — extension 
to the usual two-spin species case is straightforward and 
irrelevant for our purpose. 

The sign of Sfl can be determined by analyzing the 
relation between A(w,q) and A+((jj,q). The non-zero 
averages (TT-Cq(T)c_q(0)) and {Ttc'1^{t)c^_^{0)) can only 
be obtained if they are taken with respect to a state with 
broken U(l)-symmetry (particle number conservation). 
If an appropriate broken-symmetry mean-field Hamilto- 
nian H]\jp exists, then 

F(T,q) = (l/Z)Tr [e-'3^^'^T^e^-f^*'^Cqe~"'-^"^c_q] ,(3) 

F+(r,q) = (l/Z)Tr 



where Z — Tr[e 



-l3H^ 



A straightforward comparison 



shows that the two Green's functions are related as 

F*(T,q) =F+(T,q), or F*(w, q) = F+(-c^, q). (4) 

As a consequence of Eq. ^ the functions A((jj,q) and 
A+(ti;,q) obey identical relation and therefore the prod- 
uct A(aj, q)A+(ti>, q) in Eq. ([T]) is negative definite for the 
odd-frequency A(a;,q) producing 5D, > 0. As the result, 
one is forced to conclude that the odd-frequency super- 
conducting phase is thermodynamically unstable.— This 
conclusion can also be reached if one uses the Green's 
functions obtained in original work by Berezinskii.^ As 
a related issue, one also finds an unphysical Meissner ef- 
fect (i.e., with the negative Meissner kernel),- Later it 
was suggested that the odd-frequency state might exist 
as an inhomogeneous phase, where the order parameter 
is modulated at the microscopic levelj^i^i^ or might be a 
manifestation of some composite (even-frequency) order 
parameter /i*' 

In what follows, we show that relation ([4]) holds only 
for the ewen-frequency anomalous Green's function. Odd- 
frequency superconducting state changes this relation. 
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and Eq. (g]) is modified to F*(r,q) = -F+(T,q). The 
problem with the reasoning outhned above is that it 
assumed existence of a mean- field Hamiltonian, Hmf- 
However, to account for the retardation effects, which are 
essential in the case of the odd-frequency state, effective 
Hamiltonian language, e.g. Eq. is inappropriate. In 
other words, Hmf for odd-frequency superconductivity 
does not exist. Instead, one should consider an effective 
action that is essentially non-local in time. 

To study the superconducting phase, we represent the 
partition function of the system is a functional integral^ 



Z = / V'ipV'iP'DA*'DAe 



(5) 



with 

S = I dxiip{xi) 



V'(a;i)+/dxidx2^^^i^^ (6) 



D{xi - X2) 

dxidx2[A* {xi, X2)i){x2)ik{x-^+A{xi, x-^'ii}(x-^i){x2)\ 



where ^{x) and ip{x) are conjugate Grassmann vari- 
ables corresponding to the fermionic fields V'U'') — 
Z]q6~*'*'"Cq and V(r) = Z]qe"'''cq, with x labeling both 
spacial r and (Matsubara) time t coordinates, and ^ is 



the kinetic energy operator, ^ 



-V2/(2m) 



In 



Eqs. ([SI [6]) we have introduced the pairing field A(a;i, 0:2) 
via the standard Hubbard-Stratonovich transformatioi*ii 
by decoupling the interaction term J dxidx2D{xi — 
X2)4'{xi)i^{x2)i^{x2)'4'{xi). Note that up to this point, 
no approximation has been made. 

As pointed out earlier, the anomalous Green's func- 
tions can be defined only with respect to a state with the 
broken U{1) symmetry. Indeed it is easy to see that the 
quantity 



-S(i/',i/',A*:A) 



(7) 



is identically zero: after integration over the A fields we 
average 'ip{x)^{x') with respect to the action containing 
only products V'V' of Grassmann variables. This aver- 
age is nothing but F{x — x') defined earlier in terms of 
the time-ordered average, which is indeed zero in normal 
phase. While in the normal phase, the partition func- 
tion Z of Eq. ([5]) is dominated by the vicinity of A = 0, 
below certain temperature the situation changes: The 
primary contribution to Z comes from |A(a::i,a;2)| 7^ 0, 
which signals spontaneous U{1) symmetry breaking. In 
addition, for spinless electrons, the superconducting state 
breaks either spatial parity (e.g. p-wave superconductor) 
or time reversal symmetry (odd-frequency superconduc- 
tor). In the ordered state, one can expand the action 
in the vicinity of the non-zero saddle-point value of the 
order parameter Amf{xi — ^2), with mean field approx- 
imation corresponding to further neglecting the fluctua- 
tions around the saddle point. At the mean-field level. 



the anomalous correlation functions {F and F'^) can be 
expressed as 

F(T-r';r-r') = Z^,VyPV5W(r,r)^(r',r')e-^"-,(8) 
F+(t-t'; r-r') - Z]^/pJv^V^i,ir, r)7^(r', T')e-^''^ ,(9) 
where 

Smf = [ dxi-ipixi) dr+i iiixi) (10) 

(dxidX2A\jp{xi —X2)'4'{X2}'4'{XI) 

' dxidx2AMF{xi - X2)4>ixi)i>{x2), 



and AMF(a;i —X2) is again defined by the self-consistency 
condition, Eq. 

Now we are ready to determine the relation between 
F and . Due to the long history of the problem, and 
since we believe that this relation is the root of the di- 
vergent claims about the fate of the odd-frequency super- 
conductivity, we present here all technical details. First, 
let us take the complex conjugate of F: 



(11) 



where 



S*MF = / dx ^^dr^ix)^P{x) + i^ix)^P{x)^ (12) 



dxidx-2 



Amf{xi - X2)'4'{xi)i!{x2) 



+ A\jp{;Xi - X2)i'{x2)'4'{xi) 



Integrating the first two terms in Eq. (jl2p by parts we 

obtain J dxTp{x) ~dr + £, ijj{x). Then defining the new 

variables according to ^(r, r) = 7/'(r, — r) and ?7(r, r) = 
-(/^(r, — t), and changing r — r in every integral in 
Eq. (fT^ . we obtain 



(13) 
(14) 



F*{t -t'-v- v') = ZjlpjVf]Vr^f]{v\ -r')^(r, -T)e 
with 

Smf = / dxfi{x) dr +i r]{x) 

AA/i^(r2 - Ti,ri - r2)r](x2)f]{xi) 



dxidx2 



+Amf{t2 - ri,ri - r2)?7(a;i)77(a;2) 



For even-in-r A(r, r), (e.g., p-wave for a single spin 
species case) we have Smf = Smf and therefore com- 
paring Eqs. HSlIlll) with Eqs. ©[TOl) we recover Eq. (g]). 
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For odd-in-r A(t, r), (e.g., s-wave with a single spin 
species) we see that by changing t —t a minus sign is 
generated in the last two terms in the RHS of Eq. (fT4|) 
as compared to Eq. p^ . This difference can be re- 
moved by another change of variables 77(r, t) — > iri{r, t), 
fj(r,T) — > — j77(r,T), which is a simple gauge transfor- 
mation. We obtain again Smf = Smf- However, an 
additional factor (-1) now appears before the entire path 
integrals owning to the fact that the quantity fjfj trans- 
forms into {—fjfi) as a result of the last gauge transfor- 
mation. Therefore comparing Eqs. [H]) with Eqs. ([HI 
nop we finally obtain that in the odd-frequency case we 
have the relation 

F*{t - r', r - r') = -F+{t - t', r' - r) (15) 



or 



F*iu;,q)^F+iu,q), 



(16) 



Obviously Eq. (|T6l) holds for the even-frequency case as 
well — in which case F^{—uj,q) in Eq. ^ can be re- 
placed by F^{uj,q) for even- frequency F. As a result, 
contrary to the conclusion of Refs. 3,4;j|@, the product 
A(aj, q)A+(ti;, q) in Eq. ([T]) is positive definite both for 
odd- and even- frequency pairing and therefore 6fl < 0. 
The same conclusion can be reached by directly analyz- 
ing the mean-field free energy, given by Eqs. E]), as 
the system undergoes the phase transition. Thus we con- 
clude that the odd-frequency superconducting phase has 
free energy lower than that of the normal phase. The 
magnitude of the order parameter has to be determined 
from the self-consistency condition and is non-zero below 
the superconducting transition temperature. 

It can be also verified that due to relation pB]) odd- 
frequency superconducting phase has a positive Meissner 
kernel and therefore a physically meaningful Meissner ef- 
fect. The supercurrent and the vector potential are re- 
lated as j(k) = — A;;(k)A(k), where the (Meissner) kernel 
K, is expressed^ as 



/C(k) 



(17) 



m m?l3 

+ F{u,^p+)F+{u,^p^)\ 



Here p± = p ± k/2. The Green's functions can be easily 
obtained from Eqs. ([51 [H [TUl) . We have 



F(w,q) = 



2A7\/i?(cj,q) 



e(q)2 + 4|AMF(w,q)| 



(18) 



F+(c^,q) 



2AW(^,q) 



e(q)2+4|AMF(c^,q) 



-e(q) 



^' + e(q)'+4|AMF(cc',q)| 



(19) 



(20) 



Note that the form of Eqs. (flSl [T9| is consistent with 
Eq. (HI]), not with Eq. (g]). Had we used Eq. (g]), we 
would have obtained A^,j^(— a;,q) in Eq. p9)) as well as 
-|A(w,q)|2 in the denominators of Eqs. (fT^fTOll^ . 

As usual Eq. (|17p is divergent and we regularize it by 
subtracting /C(k) for Amf(w, q) = (obviously /C(k) = 
in normal phase) In the long wave-length limit and for 
Aj\/i?(w, q) independent of q (i.e., for pairing in s — wave 
channel) we obtain 



4|AMF(t^)|- , . 

(^2+4|Amf(^)|2)3/2- ^ ^ 



This equation obviously is positive definite. Note that 
if we had used Eq. ([4]) (which is invalid as we argue 
above), we would have obtained A+(a', q)A(aj, q) = 
A*(— tj, q)A((jj, q) in the numerator in the RHS of 
Eq. (|2ip and thus negative Meissner kernel for the odd- 
frequency case. 

While possessing similar electromagnetic properties 
(i.e., the Meissner effect) to its familiar even- frequency 
counterpart, the odd- frequency superconductor is ex- 
pected to differ from the even-frequency one in several 
important aspects. Here we mention just some of them 
qualitatively. The equal spin pairing considered above, 
leads to a gapless superconductor, with an isotropic (s- 
wave) electronic spectral function. The only other known 
example where this can happen is the s-wave supercon- 
ductor with a relatively high concentration of magnetic 
impurities; however, for odd-frequency superconductor 
this would occurs even in the clean case. The odd- 
frequency anisotropic s-wave superconductor can readily 
exceed the Pauli paramagnetic limit, and can show very 
little change in magnetic susceptibility (Knight shift) 
across the superconducting transition. These and other 
properties, in particular the manifestations of the time- 
reversal symmetry braking in such superconductors, are 
attractive direction for future detailed theoretical inves- 
tigations. 

We acknowledge valuable and stimulating discussions 
with E. Abrahams and A. Balatsky. This work is sup- 
ported by the US DOE. 
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